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The mass of the triply heavy baryon Qbbb is calculated in lattice QCD with 2 + 1 flavors of light 
sea quarks. The b quark is implemented with improved lattice NRQCD. Gauge field ensembles from 
both the RBC/UKQCD and MILC collaborations with lattice spacings in the range from 0.08 fm 
to 0.12 fm are used. The final result for the Qbbb mass, which includes an electrostatic correction, is 
14.371 ± 0.004 s tat ± 0.011 syst ± 0.001 cxp GeV. The hyperfine splitting between the physical J = 3/2 
state and a fictitious J = 1/2 state is also calculated. 

PACS numbers: 12.38. Gc, 14.20.Mr 



I. INTRODUCTION 



Hadrons containing only heavy valence quarks play a key role in the study of QCD due to their approximately 
nonrelativistic nature and their clean spectrum of narrow states. Many precise experimental results are available 
for heavy quarkonium mesons, and their analysis has contributed significantly to the understanding of the quark- 
antiquark forces OH]. However, the nature of the QCD gauge symmetry is exhibited in the most direct way in the 
baryons, where the number of valence quarks is equal to the number of colors. It is therefore of great interest to 
explore triply-heavy baryons like the £lbbb- Masses and other properties of triply- heavy baryons have been calculated 
using various quark models [3TU8] and sum rules |19j . Most potential-model calculations include only two-body forces, 
but QCD also leads to genuine three-body forces. A three-body linear confinement term was found in the static three- 
quark potential calculated nonperturbatively from lattice QCD [2D]. In perturbative QCD, three-body forces first 
arise at next-to-next-to-leading order |21) . 

Triply-heavy baryons have not yet been observed experimentally (see [5j I22H25] for experimental aspects and 
calculations of production rates in colliders), and the predictions of their properties can not yet be compared to 
the real world. The model-dependent calculations can however be tested by comparing them to nonperturbative 
first-principles calculations in lattice QCD. 

While there are several recent lattice QCD calculations of singly- and doubly-heavy baryon masses [26-31], there 
appears to be only one lattice result for a triply-heavy baryon mass in the literature so far: the Vt ccc in Ref. [32] . 
The calculation in [32] was performed in the quenched approximation, i.e. neglecting the effects of dynamical light 
quarks, and the result for Mn cec obtained there may also have significant discretization errors [33| due to the use of 
a relativistic action for the charm quark at the rather large value of am c = 0.8 (where m c is the bare charm quark 
mass and a is the lattice spacing). 

In the following, a precise calculation of the mass of the triply-bottom baryon ttbbb in lattice QCD with 2 + 1 flavors 
of dynamical light quarks is reported (a preliminary result was given in [34]). The hyperfine splitting between the 
ilbbb baryon with J = 3/2 and a fictitious J = 1/2 state composed of three heavy quarks with the same mass as the 
b quark is also calculated. The computations are performed at lattice spacings in the range from approximately 0.08 
fm to 0.12 fm. At these values of the lattice spacing, the b quark can be implemented very accurately with lattice 
nonrelativistic QCD (NRQCD) 35 -37 . Most of the calculations in this work are done on gauge field ensembles that 
were generated by the RBC and UKQCD collaborations using the Iwasaki action for the gluons and a domain-wall 
action for the sea quarks [351 EH] ■ As a test of universality, calculations are also performed on gauge field ensembles 
generated by the MILC collaboration using the Liischer-Weisz action for the gluons and a rooted staggered action for 
the sea quarks 40 . The details of the gauge field ensembles are given in Sec. |n| The mass of the flbbb is obtained by 



calculating the difference En bbb — | (E nb + 3Ex), as described in Sec. Ill This quantity has only a weak dependence 
on the b quark mass and on the light quark masses, as demonstrated by the lattice results in Sec. |IV[ In the final 



result for Mn bbb , a correction due to the electrostatic Coulomb interaction of the b quarks is included (Sec. |V[). 



II. LATTICE ACTIONS AND PARAMETERS 



The problem with heavy quarks on the lattice is that relativistic actions develop large discretization errors when 
the heavy-quark mass becomes comparable to the inverse lattice spacing. In this work, the b quarks are therefore 
implemented with lattice NRQCD [351 EH], a nonrelativistic effective field theory for heavy quarks. Unlike HQET [41] . 
which is limited to heavy-light hadrons, NRQCD can be applied to hadrons containing any number of heavy quarks. 
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The use of NRQCD to study triply-heavy baryons has already been suggested in [33] ■ The power counting in this case 
should be very similar to that for heavy quarkonia, and is given in terms of the typical heavy-quark velocity inside 
the hadron. For bottomonium, one has v 2 ~ 0.1 (in units with c = 1), and the value is expected to be comparable 
for the Ofobf, (the length scales are similar; see Eqs. (17) and (18) in Sec, |v|). 

The NRQCD action employed here includes all terms up to order v , as well as Symanzik- improvement corrections 
(which reduce discretization errors). The matching to QCD is performed at tree-level, and the action is tadpole- 
improved using the mean link in Landau gauge [421 143] . The action is equal to the one used in |37| . with 
ci = C2 = C3 = C4 = C5 = cq = 1 and C7 = cs = eg = (in the notation used there.) For comparison, some results 
obtained with the order-f 2 action, which has C5 = cq = 1 and c\ = c-i = C3 = C4 = C7 = cs = eg — 0, will also be 
given. 

The details of the gauge field ensembles used in this work are shown in Table |TJ All ensembles include the effects 
of dynamical up- down- and strange sea quarks (with m u = rrid, in the following denoted as mi). The ensembles 
generated by the RBC and UKQCD collaborations [3H [3j|] make use of a domain wall action 44 -4^] for the sea 
quarks, and the Iwasaki action (47l 08] for the gluons. The domain wall action yields an exact chiral symmetry when 
the extent L 5 of the auxiliary fifth dimension is taken to infinity. The Iwasaki action suppresses the residual chiral 
symmetry breaking at finite L5 |49j (here, L5 = 16). 

Calculations are also performed on two ensembles generated by the MILC collaboration [40], using the AsqTad 
action [50-52 in combination with the rooting procedure [40] for the sea quarks, and the tadpole-improved Luscher- 
Weisz action [53 -55 for the gluons. By comparing the results from the RBC/UKQCD and from the MILC ensembles, 
sea-quark and gluon discretization effects can be disentangled from NRQCD discretization effects [37]. 

Table [I] also shows the lattice spacings of the ensembles, determined from the T(2S) — T(15) splitting [37], and 
the values of the bare 6-quark mass in lattice units, ttrnj. The 6-quark mass is set such that the kinetic mass of the 
T]b(lS) meson agrees with experiment within the statistical errors. The main calculations in this work are performed 
directly at the physical values of arrib as determined in |37j ; additional results for other values of omj, will also be 
given to illustrate the dependence on arrib- On the coarse MILC ensemble, the value of arrib used here is 0.9% below 
the physical value obtained in [37] , 
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TABLE I: Summary of lattice parameters. The bare gauge couplings are given as P = 6/g 2 (for the RBC/UKQCD ensembles) 
and P = 10 /g 2 (for the MILC ensembles). For the RBC/UKQCD ensembles, the pion masses in lattice units were taken from 
[381 1391 156] and converted to physical units using the lattice spacings given here. For the MILC ensembles, there are taste 
splittings between the different pions [40], and the root- mean-square masses taken from [57) [581 are given. 



III. CONSTRUCTION OF TWO-POINT FUNCTIONS AND ANALYSIS 



The ground-state energies of the 77b, the T, and the Qbbb are extracted from fits to the Euclidean time-dependence of 
suitable two-point functions at zero momentum. On a given lattice gauge field configuration, the two-point functions 
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for the 77b, the T, and the £lbbb are defined as 

C^{t,t',x>) = ^G$(t,x,l/ ) a>)Gp a (t > x ) l/,x'), (1) 

X 

CW(t,t',x') = ^( 7j75 )^G^(t,x,t',x')( l5lj ) lS G b s a a (t,x,t',x'), (2) 

X 

C$l 5 (t,t',x') = ^ eabce/5 ,(C7 J 0^(^)p.G^(^, a; ,^^ a; 0G^; 9 p (^, a; ,^^a : 0GS(^, a; ,^^ : r'). (3) 



Here, a, 6, c, ... are color indices (running from 1 to 3), a, (3, 7, ... are spinor indices (running from 1 to 4), C — 7072 is 
the (Euclidean) charge conjugation matrix, and the ovcrlinc denotes the Dirac conjugate. The indices j, k are in the 
range from 1 to 3. In the nonrelativistic gamma-matrix basis used here, the NRQCD heavy-quark propagator G has 
vanishing lower spinor components (for t > t'): 



G(t,x,l/,x') = t G V> 



(4) 



where G^ is the NRQCD propagator with two spinor components. The propagator G<p may include at the source 
and/or sink a Gaussian smearing operator 



ns J 



(5) 



where A^ 2 ) is a covariant lattice Laplacian. The smearing is intended to improve the overlap with the ground state 
and reduce the contamination from excited states. For the meson two-point functions, smearing is only performed on 
G^j, not G\l, while for the baryon two-point functions all three G^'s are treated equally. 

When defined through Eq. (fol, the two-point function G^ aS (t, t', x') couples to both the physical spin-3/2 state 
tlbbb, and an unphysical spin-1/2 state (when all three quark flavors are equal, this state violates the Pauli exclusion 
principle). At large Euclidean time separation t — t', the two-point function (after averaging over gauge configurations) 
approaches the form 



Cf Z% /2 e- E ^ (*-0 i(l + 7 o)(^ - friTfc) + Zf /2 e-^> <«-0 1(1 +70)^ 



(6) 



where E3/2 = En bbb and £1/2 are the ground-state energies of the J = | and J = \ states, respectively (see for 
example [59 ). The J = | and J — \ contributions can be disentangled by multiplying with the projectors 



p!- /2) = (*«-hi7j), 



13 

j(l/2) 



which gives 



(7) 
(8) 

(9) 



Before the fitting, an average over all non-vanishing spinor- and Lorentz components is calculated, which is defined 

as 



/p(j) C (q)j 



E 



P (j) r m 



i, k, a, 5, 
[(i+7o)p< ) : ,) ] q ^o 



(l + 7o)^ J) 



(10) 



In order to increase statistics, the two-point functions are calculated for 32 different source locations (t' , x') spread 
evenly across the lattice on each gauge field configuration. The source locations are shifted randomly from config- 
uration to configuration. No significant autocorrelations were seen either between source locations or in molecular 
dynamics time. An example of a matrix of tlbbb two-point functions is shown in Fig. [I] The four different functions 
correspond to smeared or local sources and/or sinks. The data are fitted by a sum of exponentials using the Bayesian 
technique from (60) . The fit functions and and priors are chosen as discussed in |61j . 
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FIG. 1: Left panel: matrix of Q.bbb two-point functions C(t) = (p( 3 / 2 > C^) with local and smeared interpolating fields (the 
local-smeared and smeared- local data coincide); the lines are from a fit with 7 exponentials and t m i n = 5. Right panel: 
corresponding effective-energy plot and ground state energy. The data are from the RBC/UKQCD ensemble with L = 32, 
ami = 0.004. Lattice units are used. 

Due to the use of NRQCD, the energies extracted from fits of two-point functions contain a shift that is proportional 
to the number of heavy quarks in the hadron. This shift cancels in the energy differences 



aEn. 



(aE Vb + 3aEr) 



(11) 



and 



aEn b 



(12) 



As will be shown in Sec. IV the quantity (11), which contains the spin average of the IS bottomonium states, has 
a weaker dependence on the 6-quark mass than (12 1 and is therefore preferred for the determination of Mn bbb . The 



energy differences (11) and (12) are computed using statistical bootstrap in order to take into account correlations. 

After conversion to physical units using the lattice spacings from the T(2S) — T(1S) splittings, the data from the 
RBC/UKQCD ensembles are extrapolated linearly in m 2 to the physical pion mass. As in [37], the data from the 
L = 24 and L — 32 ensembles, which have the same physical box size, are extrapolated simultaneously, allowing for 
an arbitrary dependence on the lattice spacing a. Higher-order terms depending on both a and m n are neglected. 



IV. LATTICE RESULTS 



nth I 



and the unphysical spin splitting aEfi bbb — aE 1 i 2 are given 



The lattice results for the energy differences 
in Table [TT} On the RBC/UKQCD ensembles wrETTL ^"24, ami = 0.005 and L = 32, ami = 0.004, results for multiple 
values of ami, are listed. The quark-mass-dependence of (11) and (12) is visualized in Fig. [2] As can be seen there, 
the dependence on ami, is smaller for (11). The quark-mass-dependence of the hyperfine splitting aEfi bbb — aEi/ 2 is 
visualized in Fig. [3} As can be seen in the plot, the dependence is slightly weaker than l/famj), similarly to the IS" 
hyperfine splitting in bottomonium [37j . The data are described well by fits with the function A/ (ami,) + B. 

On the L = 24, ami = 0.005 and L = 32, ami = 0.004 ensembles, results calculated with the leading-order 
(order- v 2 ) NRQCD action are also given in Table [n] This action does not contain spin-dependent terms, so that 
aEQ bbb — aEi/2 = and aE-f — aE^ ~ 0, and consequently the quantities (11) and (12) become equal. Remarkably, 
for the splitting aE^ bbb — | (aE, lb + 3aEr) the results obtained at the same value of ami, with the order- -u 4 and with 

'J^ 1 cLCtjlOnS ^ 1 A" "»* Vvt r 1 Ace F It ii n 1 0^, T-f- ri n f\ v d FVif\F F It /"\v/~i no it 1 ii vrfn fi n n 1 n F l An F l-i r\ f-\rf-\s~\f /ii4 

aE nbbb and § {aE nb 4 
for the splitting aEn 



the order-ii 2 actions differ by less than 1%. It appears that there is a large cancellation of the order-w 4 corrections in 

;bb . On the other hand, 
2 actions differ by about 



3aEy), 



^aEy 



which is another reason to use (11 ) for the determination of Mq 

4 



the values obtained with the order-u and with the order- 
0.1. 



10%, as might be expected for v 2 

The results for aEfi bbb — | (aE Vb + 3aEr) and aEn bbb —aEi/ 2 from the order-?; 4 action at the physical values of ami, 
were then converted to physical units using the inverse lattice spacing values from Table [Tj The chiral extrapolations 
of the data from the RBC/UKQCD ensembles are visualized in Figs. [4] and [5j As can be seen there, the dependence 



5 



of the results on the pion mass is weak. This also indicates that the higher-order terms depending on both a and 
are small. The data from the L = 16 ensembles, which provide a test of finite-volume effects, were extrapolated 
independently. No significant finite- volume effects are seen, as expected for the small size of the Qbbb an d bottomonium 
ground states (Eqs. (17) and (18)). 
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TABLE II: Results in lattice units. The errors are statistical/fitting only. 
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FIG. 3: Dependence of the (unphysical) spin splitting aEn bbb — aE\n on the bare heavy-quark mass. Left panel: RBC/UKQCD 
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FIG. 5: Chiral extrapolation of the (unphysical) spin splitting En b 
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E1/2 from the RBC/UKQCD gauge field ensembles. 
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The results of the chiral extrapolations of Ea bbb — | (E Vb + 3Ey) and En bbb — -B1/2 to = 138 MeV are given 
in Table |III| Also shown in the table are the results from the MILC ensembles, and the results from the L = 24 
and L — 32 RBC/UKQCD ensembles after interpolation/extrapolation to = 460 MeV and — 416 MeV, 
respectively, so as to match the root-mean-square pion masses of the MILC ensembles. As mentioned above, on the 
coarse MILC ensemble the value of am& used here is 0.9% below the physical value obtained in [37] ■ However, this 
affects the hyperfine splitting in physical units by only 0.2 standard deviations, which is negligible. 

For the hyperfine splitting, no significant dependence on the lattice spacing is seen, and the results from the 
RBC/UKQCD and MILC ensembles at the matching pion masses are in agreement, indicating that gluon discretization 
errors for this quantity are small. However, for this splitting systematic errors of order a s ~ 20 . . . 30% and v 2 « 10% 
due to missing radiative and relativistic corrections in the NRQCD action are expected. The final result for the 
(fictitious) hyperfine splitting is then 

En bbb - E 1/2 = 24 ± 1 stat ± 6 syst MeV. (13) 

The quantity Ea bbb — | (E Vb + 3Ey) is seen to change by about 6 MeV (Lie) when going from the coarse to 
the fine lattice spacing on the RBC/UKQCD ensembles. At the coarse lattice spacing and at the matching pion 
masses, the results from the RBC/UKQCD and MILC ensembles differ by about 3 MeV (1.0<r); at the fine lattice 
spacing this difference is 4 MeV (1.2a). Given that the ratio of a 2 between the fine and coarse lattices is about 
2, it seems reasonable to assume a maximum total discretization error of 10 MeV (5%) for the result from the fine 
(L — 32) RBC/UKQCD ensemble. The systematic errors from the missing radiative and relativistic corrections in 
the NRQCD action can be estimated to be of order a s v 2 ~ 2 . . . 3% and v A ~ 1% respectively. Thus, the final result 
for Ea bbb — § {E Vb + 3Ey) (without the electrostatic correction to be discussed in Sec.|v]) is 

En bbb - I [E Vb + 3E T ) = 0.198 ± 0.003 stat ± 0.011 syst GeV. (14) 



Collaboration 


L 3 x T 


(GcV) 


En bbb ~U E v b +3E T ) (GeV) 


En bbb - E 1/2 (MeV) 


RBC/UKQCD 


16 3 x 32 


0.138 


0.214(11) 


25.6(2.6) 


RBC/UKQCD 


24 3 x 64 


0.138 


0.2044(44) 


24.2(1.2) 


RBC/UKQCD 


32 3 x 64 


0.138 


0.1984(29) 


23.91(99) 


MILC 


24 3 x 64 


0.460 


0.2052(22) 


23.69(65) 


RBC/UKQCD 


24 3 x 64 


0.460 


0.2022(22) 


23.33(63) 


MILC 


28 3 x 96 


0.416 


0.2008(24) 


23.41(72) 


RBC/UKQCD 


32 3 x 64 


0.416 


0.1966(24) 


23.24(84) 



TABLE III: Values of En bbb — § [E Vb + SEy) and Eci bbb — £1/2 in physical units. The errors are statistical/fitting only. The 
first three rows of the table give the results from the RBC/UKQCD ensembles, extrapolated to the physical pion mass. In the 
remaining rows, the results from the MILC ensembles are compared to the results from the RBC/UKQCD ensembles; there, 
the RBC/UKQCD data were interpolated/extrapolated to match the root- mean-square pion masses of the MILC ensembles. 



THE MASS OF THE tt b 



The result ( 14 ) is the pure QCD value and needs to be corrected to include the effects of electromagnetism. Since the 
quarks in both the £lbbb and bottomonium are heavy and are moving slowly, the dominant electromagnetic correction 
is due to the electrostatic Coulomb interaction. The Coulomb interaction is repulsive in the flbbb and it therefore 
increases En bbb . On the other hand, the Coulomb interaction is attractive in bottomonium and it therefore decreases 



tbbb 

E Vb and Ey. Using (T|r _1 |T) = (i]b\r \rjb), the electrostatic correction to (14) becomes 



^coulomb = + i*^V|-|T>- (15) 

47re r I 47reo r 

Here, r is defined as the distance between two b quarks (for the f^fcb) or the distance between the b and b (for the T). 
The expectation value (T|r _1 |T) is calculated numerically using the IS wave function obtained with the QQ-onia 



package [53] for the Cornell potential with parameters as in 



This gives 



1 



T|-|T) = 8.1 fin" 



(16) 



In Ref. [S], the mean-square distance of a heavy quark from the center of mass in the £l b bb has been calculated using 
a potential model, with the result (f2bt.& I '"cm I ^fofot>) = 0.021 fm 2 . The geometry of the system implies r = \/3?*CMi an d 
therefore 



V(^bbb\r 2 \n bbb ) = 0.25 fm. 
For comparison, in bottomonium the wave function from QQ-onia gives 



v/(T|r 2 |T) = 0.20 fm. 



(17) 



(18) 



Therefore, in the following the estimate (f2t,&f,|r 



in. 



bbb, 



= (0.8±0.4)(T|r- 1 |T) = 6.5 ± 3.2 fm -1 is used. This gives 



^coulomb = 5.1 ± 2.5 MeV. 



The full mass of the Q b hh is then calculated as 



Ma. = 



Ec 



- I (E rib + 3P T ) 



E { 



Coulomb 



LQCD 





3 






PDG 8 


Ey - E Vb 



(19) 



(20) 



Here, the first term is given by (14), and the last term (the bottomonium hyperfme splitting) is taken from the lattice 
calculation [57] : 



Ey - E Vb - 



E T - E,^ 

LQCD 

60.3 ± 5.5 stat ± 5.0 



IP tensor 



syst i 2.1 eX p 



PDG 

MeV. 



(21) 



In (21 ), the ratio of the IS" hyperfme splitting and the IP tensor splitting is used [37]. The experimental uncertainty 
in the IP tensor splitting [63] leads to the last error in (21 1. In Eq. (20), the mass of the T is taken from the Particle 
Data Group to be M r = 9.4603 ± 0.0003 GeV | 



The final result for the mass of the Qbbb is then 
14.371 ± 0.004 stat ± 0.011 syst 



± 0.001 oxp GeV. 



(22) 



VI. CONCLUSIONS 



The value for M^ bbb obtained here, Eq. (22), is compared to results from various continuum calculations in Table 



IV The results from the literature range from 13.28 ± 0.10 GeV, calculated using sum rules in [19], to 14.834 GeV 
calculated using a quark model in |17j . The nonperturbative dynamical lattice QCD calculation performed here, with 
a total uncertainty of only 12 MeV, is a valuable test of the continuum models. If an experimental result for MQ bbb 
ever becomes av ailable, it will in turn become a stringent test of the lattice calculation. 

The value (22 ) obtained here satisfies the baryon- meson mass inequality Mn bbb > %My derived in [6"4Tl66| . like most 
of the results given in Table IV with the exception of those from [5] and [T5] . 

While the lattice calculation of Mfi bbb can rule out many continuum models, the sole agreement for this value from 
a model-dependent calculation with several parameters is of course not sufficient to show that the assumptions of 
that model are correct. To establish the usefulness of a model, multiple predictions using only a small number of 
parameters need to be tested. It is planned to perform lattice calculations also for excited states of the fl bbb , which 
will give more insight into the three-quark forces of QCD. Excited states have been studied using a quark model in [9]. 
Another direction is the lattice calculation of masses of triply-heavy baryons containing charm quarks: ccc, ccb and 
ebb. Lattice NRQCD is not well suited for charm quarks due to their lower mass; instead, for example the Fermilab 
method [67] can be used for them. 
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Reference Ma bbb (GeV) 



Ponce [3] 14.248 
Hasenfratz et al. g] 14.30 
Bjorken [5] 14.76 ±0.18 

Tsuge et al. [6] 13.823 
Silvestre-Brac [9] 14.348 . . . 14.398 

Jia pj] 14.37 ±0.08 

Martynenko [B] 14.569 
Roberts and Pervin [17] 14.834 
Bernotas and Simonis [18] 14.276 
Zhang and Huang QU 13.28 ± 0.10 

This work 14.371 ± 0.004 stat ± 0.011 syBt ± 0.001 , 

TABLE IV: Comparison of results for the Qbbb mass. 
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